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We determine the general structure of the partition function of the (/-state Potts model in an 
external magnetic field, Z(G,q,v,w) for arbitrary q, temperature variable v, and magnetic field 
variable w, on cyclic, Mobius, and free strip graphs G of the square (sq), triangular (tri), and 
honeycomb (he) lattices with width L y and arbitrarily great length L x . For the cyclic case we prove 
that the partition function has the form Z(A, L y x L x ,q, v, w) = X^d=o ^ T r [(Tz,A,L y ,d) m ], where A 
denotes the lattice type, c' d ' are specified polynomials of degree d in q, Tz,A,L y ,d is the corresponding 
transfer matrix, and m = L x (L x /2) for A = sq, tri (he), respectively. An analogous formula is 
given for Mobius strips, while only Tz,A,L y ,d=o appears for free strips. We exhibit a method for 
calculating Tz,A.L y ,d for arbitrary L y and give illustrative examples. Explicit results for arbitrary 
L y are presented for Tz,A,L y ,d with d = L y and d = L v — 1. We find very simple formulas for the 
determinant det(Tz,A,L y ,d)- We also give results for self-dual cyclic strips of the square lattice. 



I. INTRODUCTION 

The g-state Potts model has served as a valuable system for the study of phase transitions and critical 
phenomena H-Q], and recently there has been considerable interest in its connections with mathematical graph 
theory @]-[8|. For two-dimensional lattices, additional insights into the critical behavior have been obtained 
from conformal algebra methods 0, [hJ. On a lattice, or, more generally, on a graph G, at temperature T and 
in an external magnetic field H, the original Hamiltonian formulation of this model is defined by the partition 
function 



Z=Y,e- PH (1-1) 



with the Hamiltonian 



(y> ' 

where o~i = 1, q are classical spin variables on each vertex (site) i € G, (3 = (fc^T) -1 , and (ij) denote pairs 
of adjacent vertices. The graph G = G(V, E) is defined by its vertex (site) set V and its edge (bond) set E; 
we denote the number of vertices of G as n = n(G) and the number of edges of G as e(G). Without loss of 
generality, we take G to be connected and we take the magnetic field to single out the spin value er, = 1. We 
use the notation 

K = (3J , h = 0H , y = e K , v = y-l, w = e h (1.3) 
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Thus, the physical ranges of v are v > for the Potts ferromagnet, and — 1 < v < for Potts antiferromagnet. 
Positive H gives a weighting that favors spin configurations in which spins have the value 1, while negative H 
disfavors such configurations. For positive and negative H, the physical range of w is w > 1 and < w < 1, 
respectively. 

The original definition of the Potts model, (| 1 . lj) and (|I.2p , requires q to be in the set of positive integers 
N+ . This restriction is removed for the zero-field Potts model by the Fortuin-Kasteleyn cluster representation 



11] 



Z(G,q,v) = « e(G V (G ' } (1-4) 

G'CG 

where G is an arbitrary graph, G' — (V,E') with E' C E is a spanning subgraph of G, and fc(G') denotes 
the number of (connected) components of G'. Because (|1.4|l does not contain any explicit reference to the 
spins {<x,} or summation over spin configurations, it allows one to define the zero-field Potts model partition 
function with q not necessarily restricted to the positive integers, N+. The zero- field Potts model partition 
function is equivalent to the Tutte (also called Tutte- Whitney) polynomial T(G,x,y), a function of major 



importance in mathematical graph theory [6j-[8j], [12|]-[14( defined by 



T(G,x,y) = £ , (1.5) 

G'CG 

where c(G') = e(G') + fc(G') — n(G') is the number of independent cycles on G' . The equivalence is 

Z(G >q ,v) = (x-l) k ^(y-ir^T(G,x,y) , (1.6) 

where x = 1 + (q/v), so q — (x — l)(y — 1). 

In order to treat the Potts model in a magnetic field for non-integral q, it is necessary to have a generalization 
of Eq. (|1.4p . that is, a formula for the partition function that does not make any explicit reference to the spins 
or any summation over the spin values, since these spin values are restricted to lie in N+. F. Y. Wu succeeded 
in constructing such a generalization, which, for an arbitrary graph G, expresses Z(G,q,v,w) as a sum of 
terms from spanning subgraphs G' of G 151 (see also [1, E|). Let us label each of the connected components 



of G' as G'i, i = l,...,fc(G'). Wu's result is [2j,l5j,ll5| 



[I5j(s 
s Q,E 



k(G') 

Z(G,q,v,w)= ve[G ' ] II (<Z-l + ™ n(G:) ) (1-7) 

G'CG i=l 

Clearly, this formula defines Z(G, q,v,w) in a manner such that q need not be in N+. Eq. (|1.7p also shows 
that Z(G, q,v,w) is a polynomial in the variables q, v, and w. In the limit h — > — oo (i.e., w — > 0), Eqs. (jl . 1|) 
and ()1.2p show that configurations in which any cr, = 1 make no contribution to Z, so that the model reduces 
to the zero-field case with q replaced by q — 1: 

Z(G,?,«,0) = Z(G,g-l,t;,l) (1.8) 

In this paper we present a method for calculating transfer matrices for the g-state Potts model partition 
functions Z{G,q 1 v 1 w) in an external magnetic field H, for arbitrary q and temperature variable v, on cyclic, 
Mobius and free strip graphs G of the square (sq), triangular (tri), and honeycomb (he) lattices with width L y 
vertices and with arbitrarily great length L x vertices. Since this method enables one to calculate Z{G, q, v, w) 
for arbitrarily great strip lengths, it complements calculations for L x x L y lattice patches based on enumeration 
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of states (e.g., 16]). Using our transfer matrix method, we determine the general structure of this partition 
function as a sum of powers of the eigenvalues of the transfer matrix, multiplied by certain coefficients that 
depend only on q, not on v or w. The result that we find exhibits some interesting differences with the form that 
has been established for the Potts model partition function on lattice strips in the case of zero external field, 
and we explain how our more general structure reduces to the zero-field form when the external field vanishes. 
We shall present explicit results for arbitrary L y given for Tz,A,L y ,d with d — L y and d = L y — 1, and the 
determinant det(Tz.A,L y .d)- We have calculated the full transfer matrices up to widths L y — 3 for the square, 
triangular, and honeycomb lattices and L y — 2 for the cyclic self-dual strip of the square lattice. Since the total 
dimensions of these transfer matrices increase very rapidly with strip width, it is not feasible to present many of 
the explicit results here; instead, we concentrate on general methods and results that hold for arbitrary L y . In 
Ref. [17] . besides mentioning briefly our structural results for cyclic lattice strips, we have used the Wu formula 
(II. 7[) to derive properties of Z(G,q,v,w), for arbitrary graphs G, concerning factorization, monotonicity, and 
zero-free regions. In Ref. we have also presented a generalization of the Tutte polynomial that corresponds 
to Z(G, q, v, w) and have formulated and discussed two related weighted graph coloring problems. Some earlier 
work using transfer matrices for the calculations of the zero-field Potts model partition function for arbitrary 



q and v on lattice strips of fixed width and arbitrary length is in Refs. 18|- 29|. Transfer matrix and related 
linear algebraic methods, as well as related generating function methods, have also been used to calculate a 
particular special case in zero field, namely the chromatic polynomial 13, Il4 , 1301 ] ; references to the literature 
can be found in reviews such as Refs. Q- [fi- 



ll. GENERAL STRUCTURE OF POTTS MODEL PARTITION FUNCTION ON LATTICE 

STRIPS IN A MAGNETIC FIELD 

A. Basic Method of Analysis and Structure for Cyclic and Mobius Lattice Strips 

In this section we derive the general structural form of the Potts model partition function Z(G S , q,v,w) 
in an external magnetic field H on lattice strip graphs G s . We label the lattice type as A and abbreviate 
the three respective types as sq, tri, and he. Each strip involves a longitudinal repetition of m copies of a 
particular subgraph. For the square-lattice strips, this is a column of squares. It is convenient to represent the 
strip of the triangular lattice as obtained from the corresponding strip of the square lattice via the insertion 
of diagonal edges connecting, say, the upper-left to lower-right vertices in each square. In both of these cases, 
the length is L x = m vertices. We represent the strip of the honeycomb lattice in the form of bricks oriented 
horizontally. In this case, since there are two vertices in 1-1 correspondence with each horizontal side of a 
brick, L x = 2m vertices. Summarizing for all of three lattices, the relation between the number of vertices 
and the number of repeated copies is 

T _ ( m if A = sq or tri or Gd in i\ 

Lx -\2m if A = Ac { ' 

Here Gd is the cyclic self-dual strip of the square lattice, to be discussed further below. 

For cyclic strips, the full transfer matrix Tz,a,l v > has a block structure formally specified by 

Ly 

Tz,a,l v = ^Y\Tz,a,l v a (2.2) 

where the product Y\Tz,A,L y ,d means a set of square blocks of the form Xz,A,L y ,d,j times the identity matrix. 
As indicated, each block is indexed by a non- negative integer d, which runs from to L y . We shall refer to 
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this as the degree of the block. The reason for this terminology is that for q > 5, this abstract submatrix has 
a dimension given by a certain polynomial cS d ' defined below in Eq. (|2.4|) , which is of degree d in q. We shall 
also refer to the \z.A,L v ,d in this block as lying in the degree-ei subspace of the full space in which the transfer 
matrix is defined. From Eq. (|2.2|) . it follows that the partition function of the Potts model, in an external 
magnetic field, on an L x x L y strip of the lattice A has the general structural form 

Ly 

Z(A,L y xL x ,cyc.,q,v,w) = c^Tr[(Tz,A,L y ,d) m ] 

d=0 

L y n zh (A,L y ,d) 



d=0 j 



with m given by (|2.ip . Here the eigenvalues Xz,A,L v ,d,j depend on the lattice type A, the strip width L y , and 
the variables q, v, and w, but not on the strip length, L x . The number of different Az,A,L„,d's in each subspace 
of degree d is given by nzh(L y , d), where we use the symbol Zh to indicate the nonzero field and to distinguish 
these numbers from the different numbers nz(L y ,d) for the zero-field case. The coefficients = c^(q) are 
polynomials of degree d in q defined by: 



,^ = Y / (-iy( 2d . j y q -l) d -i (2.4) 



c 



The first few of these coefficients are — 1, = q — 2, cS 2 ^ = q 2 — 5q + 5, = q 3 — 8q 2 + I9q — 13, and 
c^ 4 ) = (q — 2)(q 3 — 9q 2 + 2Aq— 17). For q > 5, the coefficients are positive integers and can be interpreted as 
multiplicities of the corresponding eigenvalues Xz,A.L y ,d, or equivalently, as the dimensions of the sub-blocks 
Tz.A,L y ,d in the transfer matrix Tz,a,l v - Starting from this range of q, one can then continue the expression 
(12. 3[) to arbitrary q. However, for some positive integer q values, the c^'s are negative, and hence cannot 
directly be interpreted as multiplicities of eigenvalues. For example, 

= -1 for q = 2 and d = 2 mod 3 (2.5) 

and 

5 (d) =-l for q = 3 and d = 2 or 3 mod 4 (2.6) 

For brevity, we usually suppress the argument in the notation, writing simply rather than c^ d '(q). The 
coefficients for the present case of the Potts model in a nonzero magnetic field are related to the cor- 



responding coefficients 19|, [31 1 for the zero-field Potts model according to the following equation (with 



arguments indicated explicitly) 



where 



* d \q) = C W( q -l) , (2.7) 



c 



c 



(d) 



• 3 



where U n (x) is the Chebyshev polynomial of the second kind. We have also constructed a general formula 
analogous to Eq. (|2.3|) for the Z(G,q,v,w) where G is a Mobius strip and a self-dual strip of the square 
lattice, extending our earlier results for the zero-field case in \&^- 41 1. 
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To derive Eq. (|2.3[) . let us first consider the subspace of degree d = in the transfer matrix for a strip with 
width L y vertices. This matrix is defined with respect to a given basis, and we shall refer to the configurations 
that comprise this basis as the basis elements. In addition to the basis elements of the transfer matrix for 
the zero-field case consisting of all the possible non-crossing partitions of L y vertices, there are additional 
basis elements where certain vertices are in the q — 1 state. The eigenvalues of the transfer matrix for a 
free strip are the same as the eigenvalues of the transfer matrix Tz,A,L y ,d=o in this degree d = subspace for 
the corresponding cyclic strip, where a set of horizontal edges connecting two adjacent sets of L y vertices do 
not occur. The dimension of this matrix, denoted as nzh{L y ,Q) 1 is the binomial transform (j^ ] of a Catalan 
number, 



L, 



n zh (L y ,0) = ^2( *)C k (2.9) 



fc=0 

where 

c «=M 2 t) (2 - io) 

is the Catalan number. (No confusion should result from our use of the same symbol C n to denote the circuit 
graph with n vertices since the meaning will be clear from context.) Parenthetically, we note that the binomial 
transform in Eq. (|2.9j) has been of interest in other combinatorial problems and appears as sequence A007317 
in Ref. [35I ; some other sequences riven here also have appeared in different mathematical contexts and have 
been similarly catalogued in Ref. [35(. We do not show the dependence on the lattice A explicitly in nzh(L y , d) 
as it is the same for cyclic strips of the square, triangular and honeycomb lattices as for the zero-field transfer 
matrices 



23J, |3l|. (Below we shall consider self-dual strips Gd of the square lattice, which have different 
dimensions nzh{GD, L y , d); for this case we shall include the Gd dependence in the notation.) The degree 
d = 1 subspace is given by all of the possible non-crossing partitions with a color assignment, out of q — 1 
states, to one vertex, with possible connections with other vertices, plus the basis elements where certain other 
vertices are in the q = 1 state. The multiplicity is given by cf- 1 ' = q — 2. This follows because there are q — 1 
possible ways of making this color assignment, but one of these has to be subtracted, since the effect of all the 
possible color assignments is equivalent to the choice of no specific color assignment, which has been taken 
into account in the level subspace. Equivalently, c^(q) = c^(q — 1). In this derivation and subsequent 
ones we assume that q is an integer > 5 to begin with, so that the multiplicities are positive-definite; we then 
analytically continue them downward to apply in the region < q < 5 where can be zero or negative. 
For the next subspace d = 2 we consider all of the non-crossing partitions with two-color assignments to two 
separated vertices (with possible connections with other vertices) , plus the basis elements where certain other 
vertices are in the q = 1 state. This method is then continued for higher d up to the maximum degree, d = L y . 
The multiplicity (q) for general d is given by Eq. (12.41) . 

To illustrate the method further, we list graphically all the possible partitions for the strips with L v = 1, 
L y = 2, and L y = 3 in Figs. [T]-[3l where white circles are the original L y vertices, each black circle corresponds 
to a specific color assignment, and the crosses are the vertices in the q — 1 state. In the following discussion, 
we will simply use the names white and black circles and crosses with the meaning understood. We denote 
the set of partitions Vl v ,cI for 1 < L y < 3 as follows. For simplicity, a single white circle is not given explicitly 
in the notation for a partition, where contiguous vertex numbers are in the same state (color), vertices with 
overline are color-assigned, and vertices with underline are in the q — 1 state. In the following set of partitions 
Vh ,d, individual partitions are separated by a semicolon. For each partition, vertices that are not in the same 
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FIG. 1: Partitions for the L y = 1 strip. 
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FIG. 2: Partitions for the L y = 2 strip. 



state are separated by a comma. 
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Vi.o ={l;l}, Vi,i = {1} (2.H) 

P 2 ,o = U;I;2;12;12} , 7> 2jl = {T;l,2;2; 1,2; 12} , 7\ 2 = {1,2} (2.12) 
^3,0 = {-?; I; 2; 3; 12; 13; 23; 123; 12; 12, 3; 13; 13, 2; 23; 1, 23; 123} , 
Vs,i = {3; 1, 3; 2, 3; 12, 3; 2; 1, 2; 3, 2; 13, 2; T; 1, 2; 1, 3; T, 23; 12, 3; T2; 12, 3; 13; 13, 2; 23; 
1,23; 1, 23; 123} , 

V 3 , 2 = {2,3;1,2,3;1,3;1,3,2;1,2;1,2,3;I2,3;I,23} , V 3 , 3 = {1,2,3} (2.13) 



From Eq. (12. 2|) . it follows that the dimension of the total transfer matrix, i.e., the total number of eigenvalues 
^z,A,L y ,d,j, counting multiplicities, is 

Ly Ly 

Cz.cyc.Ly = dim(Tz A , Ly ) = dim (Tz,A,L y ,d) = J2 c (d) n zh (L y , d) = q L « (2.14) 
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FIG. 3: Partitions for the L y = 3 strip. 



G 



which is the sum of coefficients, and is independent of the length m of the strip. We define Nzh,A,L y as the 
total number of distinct eigenvalues of Tz.a.l v , i.e. the sum of the dimensions of the submatrices Tz : A t L y ,d> 
modulo the multiplicity c^ d '; 

Ly 

N zh , Ly =J2 n zh(L y ,d) (2.15) 

Here we include an h after the Z in the notation Nzh.\j, „ to indicate the presence of a nonzero magnetic field 
and to avoid confusion with our earlier notation in Ref. for the different (smaller) total number N z ,A,L y 
for the zero-field case. 

Using these methods, we have determined the nzh(L y ,d). We find that they are given by the following 
theorem, whose proof is similar to that for the zero-field case given in . 



Theorem II. 1 The dimension of Tz,A.L y .d, nzh(L yi d), for H =/= and < d < L y is determined as follows. 
One has 

n zh (L y ,d)=0 for d > L y (2.16) 



n zh {L y ,L y ) = l (2.17) 

and 

n zh (l,0) = 2 (2.18) 
All other numbers nzh(L y ,d) are then determined by the two recursion relations 

n Zh (L y + 1,0) = 2n zh (L y , 0) + n Zh (L y , 1) (2.19) 

and 

nzh(L y + l,d) = nzh(L y ,d-l)+3nzh(L y ,d) + n Z h(L y ,d+l) for 1 < d < L y + 1 (2.20) 

Proof. Since the maximum number of colors to assign is L y for a strip with width L yi it follows that 
nzh(Ly, d) = for d > L y and nzh(L y , L y ) = 1. It is elementary that riz/i(l,0) = 2, as shown in Fig. 
[TJ The d = partitions of a width-(Lj / + 1) strip can be obtained by either adding a unconnected white circle 
or a cross to the bottom of the d = partitions of a width-L y strip, or converting the black circle of the d = 1 
partitions of a width-L y strip into a white circle. This gives Eq. (|2.19[) . which is a special case of the following 
discussion. For general < d < L y + 1, the partitions of a width-(L y + 1) strip can be obtained in one of 
the following four ways: (a) for 1 < d < L y + 1, adding a pair of connected circles, one black and one white, 
(but not connected to any other vertex) above the highest black circle of the d — 1 partitions of a width-Lj, 
strip; (b) for < d < L y , adding a unconnected white circle or a cross above the highest black circle of the 
d partitions of a width-Lj, strip; (c) for 1 < d < L y , adding a white circle above the highest black circle of 
the d partitions of a width-L y strip and connecting these two circles; and (d) for < d < L y — 1, converting 
the highest black circle of the d + 1 partitions of a width-L^ strip into a white circle. Now the lowest white 
circle of the d partitions of a width- (L y + 1) strip can either connect to a black circle with or without ((c) or 

(a) ) other connections to other white circles, or it may not connect to a black circle with or without ((d) or 

(b) ) other connections to other white circles. The case in which the lowest vertex is a cross is included in (b). 
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TABLE I: Table of numbers nzh(L y ,d) and their sums, Nzh,L y - Blank entries are zero. 
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Therefore, (a) to (d) exhaust all the possibilities. The numbers for these four categories are nzh{L y ,d — 1), 
2nzh{L y ,d), nzh(Ly,d), and nzh{L y ,d + 1), respectively, so we have 

nzh(L y + l,d) = nzh{Ly,d — l) + 3nzh(Ly,d) + nzh(Ly,d+l) forl<d<L y — l 
nzh(L y + l,L y ) = n Z h(L y ,L y - 1) + 3nzh.{L y ,L y ) 

n Z h{L y + l,L y + 1) = n zh (L y ,L y ) (2-21) 

Since nzh{L yi d) — for d > L y , these can be combined into Eq. (|2.20p . The relation for d = 0, i.e. Eq. 
(|2.19p . results from (b) and (d) only. This completes the proof. □ We remark that these results can also be 
obtained from Eq. (12.14|) with given in Eq. (|2.4|) by the argument in the proof of Theorem 4 in Ref. • 

A corollary is that 

nzh(Ly,L y -l) = 3L y -l (2.22) 

In Table [J we list the first few numbers nzh(L y ,d) and the total sums Nzh,L y - In particular, the numbers 
n zh(L y , 0) have been given in Eq. (|2.9[) . and the numbers nzh(L y , 1) are the binomial transforms of the first 
differences of the Catalan numbers. 



Two corollaries of Theorem III. 1 1 arc the following. First, 

N Zh ,L y +i = 5N zhtLy - 2n zh (L y ,0) (2.23) 

Second, Nzh,L y can be expressed as 

= E ( L :) (1) (2-24) 



3=0 



3 / \ 3 



The construction of the transfer matrix for each level (i.e., degree) d can be carried out by methods similar to 



those for the zero-field transfer matrix 29j. Using the basis elements described above (e.g. Eqs. (|2.11|) - (|2.13|) 



for 1 < L y < 3), we define JL y ,d,i,i+i as the join operator between vertices i and i + 1, i.e. 

J L y ,d,i,i+lVv = Vpi(i+i) , (2.25) 
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where both of these vertices i and i + 1 are in states different from the q — 1 state and have not already 
been assigned different colors. Vi(i + 1) denotes the partition with d color assignments obtained from V by 
connecting vertices i and i + 1 (regardless of whether they were already connected or not) . For each subspace 
d, we also define DL y ,d,i as the detach operator on vertex i such that 

J Vp\j + v-py; if i is a cross or it is connected to other vertices . . 

L y ,d,i v — s ^ _ i)( V p i _|_ V pj) if i is a white circle without connection ' \ ■ ) 

where the vertex i should not have been assigned color. V\i is the partition obtained from V by making i a 
white circle without connection, and similarly, V\i is the partition obtained from V by making i a cross. Vi is 
the same as V if i has originally no connection, and Vi is the partition obtained from V by converting i into 
a cross. Since the application of these operators does not increase the number of colors assigned, i.e. q — 1, 
the full transfer matrix in Eq. (|2.2[) has a triangular block (submatrix) form, and the block corresponding to 
a d-color assignment has a diagonal block form with &- d ' blocks. The transfer matrix Tz,h,,L y ,d for each d is 
the product of the transverse and longitudinal parts, Hz,A,L v ,d and Vz,h,L y ,d, which can be expressed as 

Ly-l 

Hz,sqX v ,d — Hz,triX v ,d = K Y[ {I + v JL y ,dd,i+l) 

i=l 

[Ly/2] l(L v -l)/2] 
Hz,hc,L y 4,l K I [ (I + V Ji s ,d,2i-l,2i) ) Hz,hc,L y ,d,2 = K JJ (/ + V Jl v ,d,2i,2i+l) 



i=l 



L 



y 



Vz,sq : L y ,d — Vz,hc,L y ,d — + DL y .d,i) 

i=l 

Ly-l 

V z ,tri,L y ,d = ]J[{vI + D Lytd ^{I + vJ Ly4 ^+i)](vI + D LyALy ) , (2.27) 

i=l 

where [v] denotes the integral part of v. Here K is the diagonal matrix with diagonal element w e , where £ is 
the number of vertices in the q = 1 state for the corresponding basis clement. We have 

Tz,sq,L y ,d — Vz,sq,Ly,dHz,sq,L y ,d , Tz,tri,L y ,d = Vz ,tri,L y ,dH Z ,tri,L y ,d 

Tz,hc,L y ,d = {Vz,hc,L y ,dHz,hc,Ly,d,2)(Vz,hc,Ly,dHz t hc,L v ,d,l) = Tz,hc,L y ,d,2Tz,hc,L y ,d,l (2.28) 

We explain various details for both cyclic and Mobius strips. Consider two adjacent sets of L y vertices 
and denote these as ii,Z2, ■■•iL y and ji,j2, ■■■3l v - For cyclic strips of the square lattice, the horizontal edges 
connecting these vertices are («i,ji), (12,^2), {iz y ,jL v )- For Mobius strips, one set of horizontal edges 
becomes («i, Zl h ), (*2, Jl b -i), (iL y ,3i)- This corresponds to exchanging the pair of basis elements that 
switch to each other when the vertices reverse in order, i.e., the set of basis elements that do not have self- 
reflection symmetry with respect to the center of the L y vertices. For example, among the partitions for the 
L y = 2 strip in Fig. [21 the second partition 1 and the third partition 2 in 7^2, must be exchanged under this 
reflection. Similarly, the pairs of partitions in 7*2,1 are (i) the first partition 1 and the third partition 2, (ii) 
the second partition 1,2 and the forth partition 1, 2. For this specific set of edges of Mobius strips, the pairs of 
columns of Vz,A,L y ,d that correspond to these pairs of partitions should be exchanged, and these matrices will 
be denoted as Vz,h,L y ,d- Equivalently, the same pairs of columns of Tz t A,L v ,d should be exchanged, and these 
matrices will be denoted as Tz,k,L y ,d — Vz,h,L y ,dHz,h,L y ,d for A = sq,tri. There are two kinds of Mobius 
strips for the honeycomb lattice. When L y is even, the number of vertices in the horizontal direction is even 
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as for the cyclic strips, i.e., L x = 2m. When L y is odd, the number of vertices in the horizontal direction is 
odd, L x = 2m — 1. Therefore, for the honeycomb lattice, we use the definition 



T. 



Vz,hc 



L v ,d,Hz.hc.L^.d,l 



for odd L„ 



TzMc,L y ,d — Vz,hc,L y ,dHzMc,L y ,d,2Vz,hc,L y ,dHzjic,L y ,d,i for even Ly (2.29) 

As was discussed in for the crossing-subgraph strips, the square of each eigenvalue of Tz,hc,L y ,d for odd 
L y is an eigenvalue of the corresponding Tz : hc,L y ,d- 

We now apply these general methods to determine the structure of the Potts model partition function in a 
magnetic field on a lattice strip with Mobius boundary conditions. In the case of zero external magnetic field, 
we previously determined the changes of coefficients when the longitudinal boundary condition is changed 
from cyclic to Mobius 29, 3l|. Here we have the same changes of coefficients for the square, triangular 
and honeycomb lattices, as follows: 



g(0) 



g(0) 



(2.30) 



1 < k < 



L 7, 



and 



g(2fe+l) 



Ly-1 



0<k< 

We thus find the following general structure for the Potts model partition function for Mobius strips: 
Z(A,L y xL x ,Mb,q,v,w) = dQTr[(T z ,A,L„o) m - 1 fz,A,L y ,a] 



[{L y -\)/2] 

E < 

d=0 

[Ly/2] 



- (d+1) Tr[(T ZtA r 2d+1 ) m - 1 f z , A!L o d+1 ] 



E 

d=l 



Tr[(Tz,A,L 2d) m Tz,A,L y ,2d] 



(2.31) 



(2.32) 



(2.33) 



For the square lattice or the honeycomb lattice with L y even, the eigenvalues of Tz t A,L y ,d are the same as 
those of Tz,A,L v ,d except for possible changes of signs. The number of eigenvalues with sign changes is equal 
to the number of column-exchanges from Tz,\,L y ,d to Tz,A,L y ,d- Denote the number of eigenvalues that are 
the same for Tz lS q.L y ,d and Tz lS q,L y ,d as nzh(sq, L y , d, +), and the number of eigenvalues with different signs 
as nzh(sq, L y ,d,—). It is clear that 



n Z h(L y , d) = nzh(sq, L y , d, +) + nzh(sq, L y ,d, -) 



(2.34) 



Define 



An Z h(sq, L y , d) = n Z h(sq, L y ,d, +) - n Z h(sq, L y ,d, -) 



(2.35) 



which gives the number of partitions that have self-reflection symmetry. For example, among the partitions 
for the L y = 2 strip in Fig. [21 the partitions /, 12 and 12 in 7^2, o> the fifth partition 12 in 7-2,1, and the 
partition 1, 2 in 7^2,2 have self-reflection symmetry. Among the partitions for the L y — 3 strip in Fig. 02 those 
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with self-reflection symmetry include (i) the first partition /, the third partition 2, the sixth partition 1_3, the 
eighth partition 123 , the eleven partition 13, the twelve partition 13, 2 and the fifteenth partition 123 in T-^o, 
(ii) the fifth partition 2, the eighth partition 2, 13, the sixteenth partition 13, the seventeenth partition 13,2, 
and the twenty first partition 123 in T^i, (iii) the third partition 1,3 and the fourth partition 1,3,2 in V-&,?.-, 
and (iv) the partition in T>% 3 . We list Anzh(sq, L y ,d) for 1 < L y < 10 in Table UTl The relations between 
Anzh{sq,L y ,d) are 

An Zh {sq,2n 7 0) = 2An zh (sq, 2n - 1, 0) - An Zh (sq, 2n - 2, 0) for 1 < n 

An zh (sq,2n+ 1,0) = 2An zh (sq, 2n,0) + An Zh (sq,2n, 1) for < n 

Anzh{sq, 2n, 2m — 1) = Anz/i(s<7, 2n, 2m) 

= Anzh(sq, 2n — 1, 2m — 1) + Anzh(sq, 2n — 1, 2m) 

— Anzh{sq, 2n — 2, 2m — 1) for 1 < m < n 

Anzh{sq,2n + l,m) = An Z h(sq, 2n, m - 1) + An Z h{sq, 2n, m) 

+An Z h(sq,2n,m + 1) for 1 < m < 2n + 1 , (2.36) 
where the formal quantity Anzh(0,d) = Sd,o is assumed. Closed- form expressions are given by 

An zh (sq,2n,2m) = Till 2j ) for < m < n 

ft, V./ / V./ • ">J 

{ E?-o (™) ( + i io) T^zt for even < m < 2n 
A„ z „(„,2„ + 1,„, = ^1 t ^ tor „ddl< m < 2 „ + l <" 7 » 

^ ^3— \j / Vj-(-(m— 1)/2/ 2_?-t-m+l — — 

The total number of these partitions for each L y , denoted as ANzh,L y , is 



5 iy / 2 for even L y 



ANzh,L y = X] An Zh (sq, L y , d) = j 

d=0 



3 x S^"- 1 )/ 2 for odd L„ 



(2.38) 



B. Structure of Z(G,q,v,w) for Free Lattice Strips 

Using the methods discussed above, we find that the Potts model partition function in a magnetic field, on 
a lattice strips with free boundary conditions, has the form 

Z(A,L y x L Xl free, q, v, w) = uJ jy H z ,A,L y ,o(Tz : A,L y ,o) L:c ~ 1 SL y for A= sq,tri 

Z(hc,L y x L x , free,q,v,w) = u[ y Hz,hc,L y: o,i(TzMc,L y: o)^ L:c ^ 1 ^ 2 HTz,hc:L y ,a : 2) S sl v , (2.39) 
where <5 is defined by 



1 for even L 
for odd L 



(2.40) 
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TABLE II: Table of Anzh(s<?, L y , d) for strips of the square lattice. Blank entries are zero. The last entry for each value of L. 
is the total number of partitions with self-reflection symmetry. 



L y \d 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


AN zh , Ly 


1 


2 


1 




















3 


2 


3 


1 


1 


















5 


3 


7 


5 


2 


1 
















15 


4 


11 


6 


6 


1 


1 














25 


5 


28 


23 


13 


8 


2 


1 












75 


6 


45 


30 


30 


9 


9 


1 


1 










125 


7 


120 


105 


69 


48 


19 


11 


2 


1 








375 


8 


195 


144 


144 


58 


58 


12 


12 


1 


1 






625 


9 


534 


483 


346 


260 


128 


82 


25 


14 


2 


1 




1875 


10 


873 


685 


685 


330 


330 


95 


95 


15 


15 


1 


1 


3125 



The element of the vector u^ y for the partition V is given by (q— l)' 73 ', where \V\ is the number of components 
for V that are not in the q = 1 state. The element of the vector sl v for the partition V is equal to unity 
if there are no connections in V, where certain vertices can be in the q = 1 state. As an example, = 
((<? — I) 2 , q— l,q—l,l,q—l) and sj = (1, 1, 1, 1, 0) for L y = 2. For a strip with free boundary conditions, only 
the d = transfer matrix for the strip with cyclic boundary conditions is needed. For the square lattice or the 
honeycomb lattice with L y even, the size of transfer matrix can be reduced to nzh(sq, L yi 0, +) due to reflection 
symmetry as for the zero-field case (27j. We list nzh(sq,L y ,d,+) and nzh(sq, L y , d, — ) for 1 < L y < 10 in 
Table Ell From Eqs. (|2~34)) and ([2~35| . these are given by 



n Z h{sq,L y ,d,+) = - n Z h{L Vl d) + An Z h{sq, L y ,d) 



n Z h{sq, L y ,d, 



l r 

2 



n Z h(L y , d) - An Z h{sq, L y , d) 



(2.41) 



As we did for the zero-field case in Ref. [31(, we define, for the case of nonzero field, the numbers of 
\z,sq.L v ,j for the Mobius strips of the square lattice with coefficients ±c' d - ) as nzh,Mb(sq, L y ,d : ±). We list 
nzh,Mb(sq,L y ,d,+) and n Z h,Mb{sq, L y , d, -) for 1 < L y < 10 in Table HVl With the Eqs. (jlOO"]) to (|232"]) . 
the relations between nzh(sq, L y , d, ±) and nzh,Mb(sq, L y , d, ±) are 

nzh,Mb(sq,L y ,0,±) = n Z h(sq, L y , 0, ±) + n Z h(sq, L y , 2, t) 



nzh,Mb(sq,L y ,k,±) = n Z h(sq, L y , 2k - 1, ±) + n Z h(sq, L y , 2k + 2, =f) 



1 



for 1 < k < 
The differences for each d are defined as 

An Z h,Mb(sq, L y , d) = n Z h,Mb(sq, L y ,d, +) - n Z h,Mb(sq, L y ,d, -) 



(2.42) 



(2.43) 



For the Mobius strip of the square lattice or the honeycomb lattice with L y even, the sign changes of the 
eigenvalues of Tz,A.L y .d can be considered as the sign changes of the coefficients. For these cases, the sum of 
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TABLE III: Table of numbers nzh( si l} L y ,d,±) for strips of the square lattice. For each L y value, the entries in the first and 
second lines are n^(sg, L y , d, +) and n,z^(sq, L y , d, — ), respectively. Blank entries are zero. The last entry for each value of L y 
is the total N zhtLy . 



(<*,+) 

L v \ (d,-) 


0,+ 

o,- 


1,+ 
1,- 


2,+ 
2,- 


3,+ 
3,- 


4,+ 
4,- 


5,+ 
5,- 


6,+ 
6,- 


7,+ 
7,- 


8,+ 

8,- 


9,+ 

9,- 


10,+ 
10,- 


Nzh,L y 


1 


2 


1 




















3 


2 


4 
1 


3 
2 


1 


















11 


3 


11 

4 


13 
8 


5 
3 


1 
















45 


4 


31 
20 


46 
40 


26 
20 


6 
5 


1 














195 


5 


108 

80 


189 
166 


124 
111 


44 
36 


8 
6 


1 












873 


6 


388 
343 


759 
729 


585 
555 


249 
240 


66 
57 


9 
8 


1 










3989 


7 


1535 
1415 


3220 
3115 


2733 
2664 


1389 
1341 


447 
428 


93 
82 


11 

9 


1 








18483 


8 


6215 
6020 


13748 
13604 


12700 
12556 


7260 
7202 


2794 
2736 


716 
704 


124 
112 


12 
11 


1 






86515 


9 


26178 
25644 


60015 
59532 


58964 
58618 


37216 
36956 


16300 
16172 


5054 
4972 


1088 
1063 


160 
146 


14 
12 


1 




408105 


10 


112032 
111159 


264365 
263680 


273575 
272890 


186450 
186120 


90972 
90642 


32398 
32303 


8440 
8345 


1555 
1540 


200 
185 


15 
14 


1 


1936881 



coefficients is given by 



"max ( J In r T 

^ sr^ a / t i\~(d\ Q i° r even L v . 

C z , sq ,L y ,Mb = }^ C Z,L y ,Mb,j = 2^ An zh: Mb{sq,L y ,d)c ( > = | q(Ly+1)/2 fof odd ^ (2.44) 



3 = 1 d=0 

where 



for even L„, 



dmax -\iLtii for odd L y (2 - 45) 
That is, 

An Z h,Mb(sq, 2L y - 1, d) = An Z h,Mb(sq, 2L y , d) = n Z h(L y , d) for < d < L y (2.46) 



In previous work we have given zero-field results for the determinants for various strip graphs G s (e.g., 24| ) 

N Z,G S ,\ 

detT z (G s )= H ■\ z .r, J''^- : (2.47) 
where cz,G 3 ,j is the multiplicity of Xz,G 3 ,j- I n the present context, these can be written, for cyclic strips, as 

Ly 

det(T z>A:Ly ) - Y[[det(T z ^ LyA )f d) (2.48) 

d=0 
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TABLE IV: Table of numbers nzh Mf>( s 9> L y , d, ±) for Mobius strips of the square lattice. For each L y value, the entries in the 
first and second lines are nzh,Mb( s 1} L y ,d,-\-) and nzh,Mb{ s 1^ L y , d, — ), respectively. Blank entries are zero. The last entry for 
each value of L y is the total Nzh,L y ■ 



(d, +) 


o,+ 


1,+ 


2,+ 


3,+ 


4,+ 


5,+ 


Nzh,L y 


Ly \ (d,-) 


o,- 


1,- 


2,- 


3,- 


4,- 


5,- 




1 


2 


1 










3 


2 


4 


3 










11 




2 


2 












3 


14 


13 


1 








45 




9 


8 












4 


51 


46 


6 








195 




46 


41 


5 










5 


219 


195 


44 


1 






873 




204 


174 


36 










6 


943 


816 


249 


9 






3989 




928 


795 


241 


8 








7 


4199 


3648 


1398 


93 


1 




18483 




4148 


3562 


1352 


82 








8 


18771 


16484 


7372 


716 


12 




86515 




18720 


16398 


7326 


705 


11 






9 


84796 


76187 


38279 


5066 


160 


1 


408105 




84608 


75832 


38044 


4986 


146 






10 


384922 


355007 


194795 


32583 


1555 


15 


1936881 




384734 


354652 


194560 


32503 


1541 


14 





and we shall extend these results to arbitrary width with a general, nonzero magnetic field below. 

III. PROPERTIES OF TRANSFER MATRICES AT SPECIAL VALUES OF PARAMETERS 

In this section we derive some properties of the transfer matrices Tz,A,L v ,d a t special values of q, v, and w. 
Some related factorizations were given in . 

A. v = 

From (|1.7p it follows that for any graph G, the Potts model partition function Z(G, q, v, w) at v — satisfies 

Z(G,q,0,w) = (q- l + w) n(G) (3.1) 
Since this holds for arbitrary values of q, in the context of the lattice strips considered here, it implies 

{T z>A ,L v ,d)v=o = for \<d<L y (3.2) 
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i.e. these are zero matrices. Secondly, restricting to cyclic strips for simplicity, and using the basic results 
n = L y L x = L y m for A = sq, tri and n = 2L y m for A = he, Eq. (|3.ip implies that 

lr[{l ZiA>cyc . iLy ) lv=o-^ q _ 1+w s )2Lvm for A = hc {6.6) 

With our explicit calculations, we find that all of the eigenvalues of the matrix Tz^\,L y .d=o f° r v = vanish 
except for one, which is equal to (q — 1 + w) Ly if A = sq,tri and (q — 1 + w) 2Ly if A = he. As will be seen, 
this is reflected in the property that det(Tz,A,L y ,d) has a nonzero power of v as a factor for L y > 2 for all 
of the lattice strips considered here. The restriction L y > 2 is made because the strips of the triangular and 
honeycomb lattice are only well-defined without degenerating for L y > 2. In the case of the square lattice, for 
the case L y — 1, the transfer matrices Tz, S q,i,o a lso has one eigenvalue q — 1 + w and the other one equal to 
zero at v — 0, and Tz )S g,i,i is the scalar v. 



B. v = -l 

The special case v — — 1 defines two new types of weighted graph coloring problems, as we have discussed 
in We recall that the chromatic polynomial P{G,q) counts the number of ways of assigning q colors 

to the vertices of a graph G such that no adjacent vertices have the same color. This "proper g-coloring" 
of the vertices of G is equivalent to Z for the zero-temperature, zero-field Potts antiferromagnet, v = — 1: 
P(G, q) = Z(G, q, — 1). Here we have a generalization of this to a weighted proper g-coloring of the vertices of 
G, as described by the polynomial 

Ph{G, q, w) = Z(G, q, -1, w) (3.4) 

For H < 0, i.e., < w < 1, this is a weighted graph coloring problem in which one carries out a proper 
(/-coloring of the vertices of G but with a penalty factor of w for each vertex assigned the color 1. For H > 0, 
this is a second type of weighted graph coloring problem, namely a proper vertex q coloring with a weighting 
that favors one color. Since this favoring of one color conflicts with the strict constraint that no two adjacent 
vertices have the same color, the range w > 1 involves competing interactions and frustration. In the limit 
w — > oo, it is impossible to satisfy the proper coloring constraint, and this is embodied in the analytic result 
that for large positive w, Z(G,q,v,w) ~ (v + iy( G )w n ( G \ which vanishes as v — * — 1. In this v = — 1 special 
case, there are reductions in the ranks of the transfer matrices Tz,A,L y ,d, i-e., some of the eigenvalues vanish. 
This yields a new set of dimensions of matrix blocks, nph(L y ,d) < rizh{L y ,d). This is a strict inequality, 
i.e., nph{L yi d) < nzh(L y ,d), for all cases except d = L y , where nph(L yi L y ) = 1 = nzh{L y ,L y ). We 
have calculated these nph(L y ,d) and have obtained a number of interesting properties of the weighted graph 
coloring polynomial Ph(L y , d). These are beyond the scope of the present work and hence will be presented 
elsewhere. 



C. g = 

By substituting q = in (|1.7j) and noting the factorization u> n ( G ») — 1 = (w— 1) X^fco^ 1 wi ' we obtain the 
result that Z(G, 0, v, w) contains a factor of (w — 1). 
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D. q = 1 



Evaluating Eq. for ? = 1, one sees that the Kronecker delta functions S aiaj — 1 for all pairs of adjacent 
vertices consequently, 

Z(G, 1, w) = e Ke ^ +hn ^ = (v + iy(G) w n(G) (3 5) 



The coefficients evaluated at q — 1 satisfy 31 1 



C^(q= l) = (-l) d (3.6) 

Hence, in terms of transfer matrices, we derive the following sum rule for the present cyclic lattice strips 
G = A,L y x L xi cyc. 

E Tr[{T z ^ Ly4 ) m ] - J2 Tr[(T z ^ Ly4 r] = (v + 1)^) W «(G) for q = 1 (3J) 

0<d<L y , d even l<d<L y , d odd 

Here the number of edges e(G) for each type of cyclic strip is 

{{2L y -l)m ifA = sg 
(3L y - 2)m if A = tri , . 

(3L„-l)m ifA = /ic 
2L y m ifA=G/j 

where m is given in terms of L x by Eq. (|2.1[) . Since it applies for arbitrary m, the sum rule (|3.7|) implies 
relations between the eigenvalues of the various transfer matrices Tz,A,L y ,d- 

E. tu = 1 

An interesting question concerns how Z(G,q,v,w) reduces when H — ► (i.e., u> — ► 1). As will be evident 
from our explicit calculations, for generic q and v, various Xz,l y ,d,j' s m a given degree-d subspace become equal 
to Xz,L y ,d',j ,s m a subspace of different degree, d! . This process gives rise to "transmigration" of Xz,L y ,d,j's; as 
we combine the term(s) <f d '[Az,L y ,d',j] m from the degree-d subspace(s) with the term c^ d ^[\z.L y ,d' .j\ m in the 
degree-d' subspace, this has the effect of yielding the zero-field term ) [Xz,L y ,d',j] m - Associated with this 



transmigration process, the nzh{L y , d) are changed to the nz{L y ,d) given in Theorem 4 and Table 3 of 31 1 



(see also [19(). 



F. w = 

The special case w = is described by the relation (|1.8[) . which is valid for an arbitrary graph G. From 
(|1.4p . it follows that Z(G, q, v) contains a factor of q. From (|1.8p it therefore follows that Z(G, q, v, 0) contains 
a factor of (q — 1). 

IV. GENERAL RESULTS FOR CYCLIC STRIPS OF THE SQUARE LATTICE 

In this section and the subsequent ones we present general results that we have obtained for the Potts model 
in an external magnetic field, valid for arbitrarily large strip width L y (as well as arbitrarily great length) for 
transfer matrices and their properties. We begin with strips of the square lattice. 
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A. Determinants 



Wc find 



det(T ZiSq , Lyid ) = (v L «) n ^ L y^ 



-i n zh {L y -l,a 



(4.1) 



where nzh(L y ,d) was given in Theorcm lII.ll This applies for all d, i.e., < d < L y with nzh(L y — l,d) = 
for d > L y — 1. The factor of u> in Eq. (|4. 1|) originates from the diagonal matrix K in Eq. (|2.27|) . and the 
power of w is the sum of the number of vertices in the q = 1 state of all the {L y , d)-partitions, which is the 
same as the power of (1 + q/v). 

Next, taking into account that the generalized multiplicity of each \z. S q,L v .d is c^ d \ we have, for the total 
determinant, 



det(T z ,sq,L y ) = X\[det(Tz,sqX v .d, 

Ly 

d=Q 



L v -1 



w^y 1 + - (l + «r* 



nzh(i H -l^)£ <d) 



(4.2) 



Using Eq. (|2.14j) together with n zh (L y , d) = for d > i y , so that ^d=o n zh{L y - 1, d)c (d) = J2d=o 1 n zh{L y 
1, d)c^ d \ we have 



drf(T WB ) = [zA]^ (l + i) Ly (1 + z;)^- 1 



(4.3) 



This determinant of the transfer matrix for the cyclic strip of the square lattice applies for arbitrary width 
L y . It is the generalization, to H ^ 0, of the zero-field result given in with the extra w factor. 



Eigenvalue for d = L y for A = sq,tri, he 



It was shown earlier that the A's for the zero-field Potts model partition function are the same for a given 
lattice strip with cyclic, as compared with Mobius, boundary conditions 
observed earlier for the special case of the chromatic polynomial, v = — 1 
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Indeed, this had been 



37j-[41| (and it was shown that 



the A's for a strip with Klein bottle boundary conditions are a subset of the A's for the same strip with torus 
boundary conditions 36|,|42j). From Theorem lII.ll one knows that there is only one A for degree d = L yi which 
we denote as \z,a,l v ,l v - That is, for this value of d, the transfer matrix reduces to 1 x 1, i.e. a scalar. We 
found that for a cyclic or Mobius strip of the square, triangular, or honeycomb lattice with width L y , 



A* 



v v for A = sq, tri 



(4.4) 



Az,A,£„,L„ = v ' 2Ly for A = hc 

which are the same as those for the zero-field case 1291. 



(4.5) 
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C. Transfer Matrix for d = L y — 1, A = sq 

From Eq. (|2.22[) it follows that for A = sq, tri or he, the number of \z,A,L y ,d,j> f° r d = L y — 1 is nzh(L y , L y — 
1) = 3Ly — l, i.e. the transfer matrix in this subspace, Tz,A,L y ,L y -i, is a (square) (3L y — l)-dimensional matrix. 
For L y = 1, 

t ( v + q-1 w \ „s 

IV,. 9 ,i,o = -i ri , x (4-6) 

y g — 1 w(l +v) J 

For L y > 2 we find the following general formula. 

(T z , sg ,L B ,L B -i)jj ■ =u L ''- 1 (27; + g- 1) for j = 1 and j = 2L y - 1 (4.7) 

(T z ,sq,L v ,L v -i)2j-i,2j-i = v L *~ 1 (3v + q-l) for L y > 3 and 2 < j < L v — 1 (4.8) 

(Tz, 3q ,L v ,L v -ihj,2j =v L *- 1 w(l + v) for l<j<L„ (4.9) 

(rz, 5 ,,V„-i)i+i,i=« irl (" + (!-l) for j= 1 and j = 2L„ - 1 (4.10) 

(T z ,^ Lb , Lh „i) 2 ,, 2 ,_i ^w^" 1 (2^ + 9-1) for L v > 3 and 2 < j < L y - 1 (4.11) 

{T z , sq ,L y ,L v -ih 3 -i,2 3 = v L y- x w for 1 < j < L y (4.12) 

{Tz,sq,L y ,L y -l)2j-l,2j+l = {Tz,sq,L y ,L y -l)2j,2j+l 
= {Tz,sqX y ,L y -l)2j+l,2j-l = (T Z ,sq,L y ,L y -l)2j+2,2j-l = V Ly for 1 < j < Ly - 1 (4.13) 

(Tz, S g,L s ,L B -i)jj = v L y (1 + «) for 2L y + 1 < j < 3L y - 1 (4.14) 

(Tz,sq,L y ,L y -lJ2j-l,2L y +j = (Tz,sq,L y ,L y -l)2j,2L y +j 
= (T Z ,sq,L y ,L v -lhj + l,2L y +j = (T Z ,sq,L y ,L y -l)2j+2,2L y +j = U Ly_1 (l + v) for 1 < j < L y - 1 

(4.15) 

(T Z , sq ,L y X y -l)2L y +j,2j-l = (T Z , sq ,L y ,L y -l)2L y +j,2j+l = V L * + 1 for 1 < j < L y - 1 (4.16) 

with all other elements equal to zero. Thus, Tz, S q,L y ,L y —i consists of four submatrices: 

1. an upper left square submatrix with indices i,j in the ranges 1 < i,j < 2L y and nonzero elements given 

by Eqs. g3D-g33i) 

2. a lower right square submatrix with indices in the ranges 2L y + 1 < i, j < 3L y — 1 and nonzero elements 
given by Eqs. (14.141) 
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3. an upper right rectangular submatrix with nonzero elements given by Eq. (|4.15[) 

4. a lower left rectangular submatrix with nonzero elements given by Eq. (|4. 16[) 

For general q, v and w, Tz, S q,L y ,L„-i has rank equal to its dimension, 3L y — 1. We illustrate these general 
formulas for the cases L y = 2 and L y = 3. For this purpose we introduce the abbreviations 



j + v 



jv - 



(4.17) 



where j is a positive integer. We have 



and 



Tz,sq,'. 



( X2 W V V\ \ 

X\ WV\ V Vl 

V X2 W Vl 

V Xl WVi Vl 

\v 2 v 2 vvi ) 



(4.18) 



T z , 



s<j,3,2 







w 


V 











''i 









Xi 


WVi 


V 











''i 









V 





X3 


w 


V 





''i 


Vl 






V 





X2 


WVi 


V 





Vl 


Vl 












V 





X2 


w 





Vl 












V 





Xl 


WVi 





Vl 






V 2 





V 2 











VVl 







V 








v 2 





v 2 








VVl 


) 



(4.19) 



In general, neglecting the v v factor, the upper left-hand submatrix has a main 2x2 block diagonal 



with end blocks equal to 



2v + q — 1 w 
v + q — 1 w(l + v) 



and interior blocks equal to 



Sv + q-1 



Adjacent to this main block diagonal are two block diagonals equal to 




2v + q-l w(l +v) J 
and the rest of the submatrix 



is comprised of triangular regions filled with O's. The upper right-hand submatrix has a band of two 2x1 
/ l + w X 



diagonals 



1 



together with triangular regions filled with O's. The lower left-hand submatrix has a band 



of two 1x2 diagonals ( v 2 ^ together with triangular regions filled with O's. And finally, in the right-hand 
lower submatrix the entries on the main diagonal are equal to v(l + v) and the rest of this submatrix is made 
up of triangular regions of O's. For the lowest values L y = 1,2, some of these parts, such as the triangular 
regions of zeros, are not present. 

Having determined the general form of Tz,sq,L v ,d for d = L y — 1, we find that a pair of its eigenvalues 
^z,sq,L y ,L y -i ,j are roots of the following quadratic equation, 



x 2 -v Ly 1 (vw + w + v + q - l)x + wv 2Ly 1 (v + q) = 
As corollaries of our general result for Tz,sq,L v ,L y —i w e calculate the trace and determinant. 

det(T z , sq ,L v ,L y -i) = uMs^-D^ (i + + v f y -i 

which is a special case of Eq. (|4.ip . and 

Tr{T z , sq ,L y ,L y -i) = v^- 1 \(q - l)L y + v 2 {L y - 1) + v{4L y - 3) + w(l + v)L y 



(4.20) 

(4.21) 
(4.22) 
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V. GENERAL RESULTS FOR CYCLIC STRIPS OF THE TRIANGULAR LATTICE 

A. Determinants 



We find 



det{T z ,tri, Lll ,d) = (v L y) n ^i L y^ 



nzh(L y —l,d) 



(5.1) 



Comparing Vz, S q,h y ,d and Vz,tri,L y ,d m Eq. (|2.27p . one sees that a set of (/ + v J l d,i,i+i) has been included 
for the triangular lattice, so that the power of (1 + v) becomes twice of the corresponding power for the square 
lattice. 

Taking into account that the multiplicity of each Xz,tri,L y ,d,j is c^ d \ it follows that the total determinant is 



det(T z , tritLy ) = f[[det(T z ^ Lytd )f d) = [v^ \w L - (l + + vf^' 1 



d=0 



(5.2) 



B. Transfer Matrix for d = L y — 1, A = tri 



From Eq. (|2.22[) . we know that the dimension of this transfer matrix is 3L y — 1 as before. The first nontrivial 
case is L y = 2. For L y > 2 we find the following general formula. 



Z,tri.L v ,L V — 1 



)i,i = v Ly -\v z + Av + q-l) 



(5.3) 



(T. 



Z,tri,L y ,L y — 1 



2j-i,2j-i =v L y 1 (v 2 + 5v + q-l) for 2 < j < L v - 1 



(T Z ,tri,L y ,L y -l)2L y -l,2L y -l = V » (2v + 



(T z ,t n ,L v ,L y -ihj,2j = v Ly : w(l + v) for 1 < j < L y 



Z,tri.L v ,L V — 1 



) 2 ,i =v L y~ 1 {v 2 + 3v + q-l) 



(T. 



Z.tri.Lu ,Lu — 1 



) 2 j,2j-i =v h »- 1 (v* +Av + q-l) for 2 < j < L y - 1 



(5.4) 
(5.5) 
(5.6) 
(5.7) 
(5.8) 



(T Z ,tri,L y ,L y -l)2L y ,2L y -l = V Ly 1 (v + q - 1) 



(5.9) 



(Tz,tri,L y ,L y -l)2j+l,2j-l — (Tz,tri, L y ,L y -1 )2j+2,2j- 1 — « Ly for 1 < j < Ly — 1 



(5.10) 



(Tz,tri,L u ,L y -i)j,k = v Ly : u> for fc = 2,4, 2£ v and 1 < j < k - 1 



(5.11) 



(Tz,tri,Ly,L y -i)j,k = v h y~ l (2v 2 +6v + q- 1) for fc = 3,5,...,2i 1/ -3 and l<j<fc-l (5.12) 



(Tz,tri,L„z y -i)j,a£„-i =v L v- 1 (v 2 +3v + q-l) for 1 < j < 2L y - 2 



(5.13) 
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{T Z ,tri,L y ,Ly-lhj+l,2Ly+j = {T Z ,tri,Ly,Ly-l)2j+2,2L y +j = V Ly 1 (l+v) for l<j<L y -l (5.14) 



{T z ,tri,L y ,L v -i) j ,2L y +k = v L *- 1 {v + l){v + 2) for 1 < k < L y - 1 and 1 < j < 2k 



(5.15) 



{T Z ,tri,Ly,L y -l)2L y +j,2j-l = V L « + 1 (3 + V) for 1 < j < Ly - 1 



(5.16) 



{T z ,tri,L y ,L y -i)2L y +j,2k-i = v L « (2v 2 + 6u + q - 1) for 2 < fc < L„ - 1 and 1 < j < k - 1 (5.17) 



(Tz,tri,Ly,L y -l)2Ly+j,2Ly-l = V L y(v 2 +3v + q-l) for 1 < j < Ly - 1 



(5.18) 



(TX t „, LyiLy _i) 2 L H+ j,fc =« Lh <« for fe = 4, 6, 2L V and 1 < j < k/2 - 1 



(5.19) 



(T Z ,tri,Ly,Ly-l)2Ly+j,2L y +k =V L «{V + l)(v + 2) for 1 < j < k < L y - 1 



(5.20) 



with all other elements equal to zero. Thus, Tz.tri,L y ,L y ~i can again be usefully viewed as consisting of various 
submatrices. 

For general q, v and w, T z ,tri,L y ,L y -i nas rank equal to its dimension, 3L y — 1. We illustrate these general 
formulas with some explicit examples for L y = 2 and L y — 3. For compactness of notation, we use the 
abbreviations 



yj = v 2 + jv + q - 1 , Zj = 2v 2 + jv + q — 1 
where j is a positive integer. Then 

( Hi W ?/3 W V1V2 \ 

V3 wvi 2/3 w VlV 2 



T Z ,tri,2,l — V 



V X2 W V\ 

V Xi WV\ Vi 

\ v 2 V3 vys wv W1V2 / 



(5.21) 



(5.22) 



Z,tri,3,2 



( ?/4 W Zq W ?/3 W V1V2 V1V2 \ 

2/3 wvi z e w 2/3 W V1V2 V1V2 

v 2/5 w 2/3 W Vi V1V2 

V 2/4 wv l V3 w v l Vl1>2 

v X2 w 

t> xi wvi ui 

W 2 ^3 VZq WV UJ/3 WW t , Wlt l 2 VV1V2 

\ w 2 «3 1*2/3 ffif2 / 

As corollaries of our general result for T z .tri,L y ,L y -i we calculate the trace and determinant. 



v 



L y (3L y — l) 



w ^[l + ^ L \ 1+v f(L y 



-1) 



det{T. 

which is a special case of Eq. (|5.1[) . and 

Tr(Tz,t rilLy ,L y -i) = v^- 1 \(q - l)L y + v 3 (L y - 1) + Av 2 {L y - 1) + v{7L y - 6) + w(l + v)L y 



(5.23) 



(5.24) 



(5.25) 
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VI. GENERAL RESULTS FOR CYCLIC STRIPS OF THE HONEYCOMB LATTICE 



Determinants 



We find 



\L V -1 



nzh.(L y -l,d) 



(6.1) 



This can be understood as follows: by an argument similar to that given before, the power of (1 + v) is 
the same as for the square lattice case. Comparing Tz yS q,L y ,d and Tz,hc,L y ,d in Eq. (|2.28[) . one sees that 
Vz,hc.L y .d — Vz, S q.L y ,d has been multiplied twice for the honeycomb lattice, so that the powers of v and 
(1 + q/v) become twice of the corresponding powers for the square lattice. In Eq. (|2.27p . both Hz,hc,L y ,d,i 
and Hz.hc,L y ,d.2 include the matrix K so that the powers of w also become twice of the corresponding powers 
for the square lattice. 

Taking into account that the multiplicity of each Xz,hc,L y ,d,j is c^ d \ it follows that the total determinant for 
the he lattice is 



det{Tz M .L v )^\[[det{Tz M .L v ,dt d) = (v 2L ^f v \w^(l + + v) L ^ 



d=0 



(6.2) 



Summarizing the connections between the determinants of the transfer matrices for the three lattice strips, 
det(T z ,tri,L y ,d) is related to det(T ZtSq x y ,d) by the replacement of (1 + v) by (1 + v) 2 and det(T Zyh c,L y ,d) is 
related to det(Tz, S q.L y ,d) by the replacements of the respective factors wbyw 2 ,i)by v 2 and (1 + -) by (1 + ^) 2 . 
This, together with the fact that nzh(L y ,d) is the same for all of these three lattices means that the total 
determinants det{Tz t tn.L y ) and det{Tz t hc.L y ) are related to det(Tz, S q,L y ) by the same respective replacements. 

B. Transfer Matrix for d = L y — 1, A = he 



From Eq. (|2.22p . we know that the dimension of this transfer matrix is again 3L y — 1. The first nontrivial 
case is L y = 2. Recall in Eq. (|2.28[) the transfer matrix for the honeycomb lattice, Tz,hc,L y ,L y -i, is the product 
of T z ,hc,L y ,L y -i,i and T z ,hc,L y ,L y -i,2- For L y > 2 we find the following general formula. 



(T z ,h- 



C 7 Ly jLy ljl 



)2i-i,2i-i = V L * 1 (2v + q-l) for 1 < j < L y - 1 



(6.3) 



(Tz,hc,Ly,Ly—l,lJ2Ly-l,2Ly-l — 



v Ly 1 (2v + q-l) for even 
v Ly ^ 1 (v + q — 1) for L y odd 



(Tz,h, 



c,L v .L y — 1.1 



2j,2j = v Ly 1 w(l + v) for 1 < j < L y 



(6.4) 
(6.5) 



Z,hc,L y ,L y — 1,1 



hj,2j-i =v L »- 1 (v + q- 1) for l<j<L y -l 



(Tz,hc 



Ly -.Ly 1.1 ) "2L y .^Ly 1 



__ J v y (v + q — 1) for L y even 



v Ly 1 (g — 1) for L y odd 



{Tz,hc,L v ,L v -i,i)2j-i,2j — v Ly 1 w for 1 < j < Ly 



(6.6) 
(6.7) 
(6.8) 
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{Tz,hc,L y ,L y -l,l)4,j-l,4,j-3 — {Tz,hc,L y ,L v -l,l)4j,4j-3 

= {T z ,he,L v ,L v -i,i)4j-a,4j-i = {Tz,hc,Ly,Ly-i,i)4j-2,4j-i = v L * for 1 < j < [L y /2] (6.9) 

(^,^^-1,1)2^+2,-1,2^+2,-1 = v L y (1 + v) for 1 < j < [L y /2} (6.10) 

{Tz, hc ,L v ,L u -l,lhL u+ 2 j ,2L u+ 2 J =V L y for 1 < j < [(L y - l)/2] (6.11) 

{Tz,hc,L v ,L v -l,l)4,j-3,2L v +2j-l = (7z ! ft C; L !; ,Lj / -l,l)4j-2,2L !; +2j-l 
= (Tz,hc,Ly,L y -l,l)4,j-l,2L y +2j-l = {Tz,hc,L y ,L y -l,l)4j,2L v +2j-l 

= v L y~ l {v + l) for l<j<[L y /2] (6.12) 

{Tz,hc,Ly,L v -l,l)4j-l,2Ly+2j = {Tz,hc,Ly,L v -l,l)4j,2Ly+2j 
= {Tz,hc,Ly,Ly-l,l)4j+l,2L y +2j = {Tz,hc,L v ,Ly-l,l)4j+2,2L y +2j 

= v^- 1 for 1 < j < [(L y - l)/2] (6.13) 

(Tz i / lc ,L I) ,i, J ,_i i i)2i, J ,+2j-i,4j-3 = (^,^^,,^-1,1)2^+2,-1,4,-1 = for 1 < j < [L y /2} (6.14) 

(Tw„i,-i,2)i,i = w^"^ + 9-1) (6.15) 

(T z , hc ,L B ,L s -i,2)2j-i,2j-i =w L ^ 1 (2w + g-l) for 2<j<L„-l (6.16) 

^ / w L « _1 (2u + q - 1) for L y odd f«i7\ 

(^0,^-1,2)2^-1,2^-1 - j vL y. 1{v + g _ 1} for £ J eyen (6.17) 

{T z ,hc,L y ,L v -i,2)2j,2j = v^wiv + 1) for \<3<L y (6.18) 

(^0,^-1,2)2,1 = w^-^g - 1) (6.19) 

(Tz,hc,Ly,Ly-i,2)2j,2 j -i=v L v- 1 {v + q-l) for 2<j<L y -l (6.20) 

rT , s / v L y^ 1 (v + q - 1) for L y odd / fi01l 

:/ '-'.'. ' -^'. -^ ' - ( v i»-i(g - 1) for L y even (6 ' 21) 

(^,^0,^,^-1,2)2,-1,2, = lA" 1 ™ for 1 < 3 < L y (6.22) 

{Tz,hc,L y ,Ly-l,2)ij+l,4,j-l = (Tz,hc,L v ,L y -ia)4j+2Aj-l 
= (Tz,hc,L v ,L y -l,2)4j-l,4j + l = {1 !z, he, L y , L y -1,2) 4 ,', 4 ,+1 

= w L « for 1 < j < [(Lj, - l)/2] (6.23) 

(T Z j lc ,L !( ,L y -i,2)2L y +2,-i,2L y +2,-i = *> L " for 1 < j < [L y /2] (6.24) 
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(Tz,fccL„L B -l >2 )2 is +2J,2L,+« = v L * (1 + v) for 1 < j < [(L y - l)/2] 



(6.25) 



(T'z,ft,c,L H ,L y -l,2)4j-3,2L y +2j-l — (7 1 Z,hc,L y ,L H -l,2)4j"-2,2L y +2j-l 
{Tz,hc,Ly,L v -l,2)4,j-l,2L v +2j-l = (Tz,hc,Ly,Ly-l,2)4j,2Ly+2j-l 

v L y- x for l<j<[£ a /2] 



(6.26) 



(Tz,hc,Ly,L v -l,2)4j-l,2Ly+2j — (Tz,hc,Ly,Ly-l,2)4j,2Ly+2j 
(Tz,hc,Ly,L y -l,2)4j+l,2L y +2j = (Tz,hc,L v ,Ly-l,2)4j+2,2Ly+2j 



(Tz, 



he, Ly , L y 



1,2) 



2L„+2j,4j-l 



(1 + w) for 1 < j < [(L y - l)/2] 



(Tz,hc,L y , L y -l,l)2L v +2j,4j+l — V L 



(6.27) 

for 1 < j < [(Ly - l)/2] (6.28) 



with all other elements equal to zero. 

For general q, v and w, Tz,hc,L v ,L v -\ has rank equal to its dimension, 3L y — 1. We illustrate these general 
formulas for the cases L y = 2 and L y = 3 (with the same abbreviations as before): 



( X\ w 
q — 1 uwi 



'-Z,hc,2,l,l 



Z,sq,2,\ 



Z,hc,3,2,l 



T z . 
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(6.29) 



(6.30) 



The Tzjic,L y ,L y - 



are obtained via Eq. (|2.28|) from these auxiliary matrices. For example, 



I xi(3v + w) + (q — l) 2 w(x\ + wv\) 

11(11 + ami) w(wvf + q 

v(x2 + w) 

v(xi + wvi) 



V 







v(x2+w) vi(x 2 + w) \ 

1) «(zi+«)Di) Hi (n — 

xi(3u + ui) + (g — l) 2 w(x\ + ioui) vi(x2 + w) 
2:1(^1 + itfUl) «i(«it) 2 + g — 1) di(xi + wui) 

) 
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(6.31) 

We find that in general, for degree d = L y — 1, a pair of its eigenvalues Xz,hc,L v ,L y -i,j are roots of the 
following quadratic equation, 

x 2 - v 2L »- 2 (w 2 (v + l) 2 + 2(q -l)(v + w) + (q - l) 2 + v 2 ^x + w 2 v AL y- 2 (v + q) 2 = (6.32) 

The expressions for the other 3L y — 3 eigenvalues are, in general, more complicated. 

VII. GENERAL RESULTS FOR CYCLIC SELF-DUAL SQUARE-LATTICE STRIPS 

In this section we consider the Potts model for families of self-dual strip graphs of the square lattice with 
fixed width L y and arbitrarily great length L x , having periodic longitudinal boundary conditions, such that 
all vertices on one side of the strip, which we take to be the upper side, are joined by edges to a single external 
vertex. A strip graph of this type will be denoted generically as Gd and, in more detail, as Gd{L v x L x ). The 
family of Gd graphs is planar and self-dual. We recall that for a planar graph G p i, one defines the (planar) 
dual graph G* ; as the graph obtained by replacing each vertex (face) of G p i by a face (vertex) of G* t and 
connecting the vertices of the resultant G* ; by edges when the corresponding faces of G p i have a common 
edge. The graph G p i is self-dual if and only if G p i = G* t . For zero-field, it is known that 



Z{G pU q, v,w = l) = v< G ^q-< G ^Z(G* pl ,q, ?-,w = l) (7.1) 



for a planar graph G p i. In general, the graph Gd(L v x L x ) has n = \V\ — L x L y + 1 vertices, equal to the 
number of faces, /. One motivation for considering the Gd strip graphs is that they exhibit, for any L y , the 
self-duality property of the infinite square lattice so that the zero-field partition function is invariant under 
v — > q/v b y (17.11) . aside from a prefactor. 

In Ref. [321 ] we gave the general form for the zero-field Potts model partition function Z(Gd, L y xL Xl q,v,w = 
1). In our current notation with m given in terms of L x by Eq. (|2.ip . this is 

Ly+1 

Z(G D , L y x L x , q, v, 1) = K (d) Tr{(T ZiGD , Lyid ) m } (7.2) 

d=l 

where 

V«U 2d -M)= d t(-iy( 2d -'- j y-i. (7.3) 
j=0 \ 3 / 



To construct the transfer matrix Tz.g d .l ,d for each d, we begin with partitions with L y + 1 vertices, 
where the single external vertex is considered as the {L y + l)-th vertex, so that the size of the matrix is 
nzh{L y + l,d) with < d < L y + 1. When the field is non-zero, the transverse and longitudinal parts, 
Hz,G D ,L y ,d and Vz,G D ,L y ,d, of the transfer matrix T z ,G D ,L v ,d can be expressed as 

Ly Ly 

Hz,G D ,L v ,d = KY[(I + vJ Ly ,d,i,i+l) , Bz,Gn,L y ,d= ^YLi 1 + vJ L y ,d,i,i+l) 
i=l i=l 

V Z ,G D ,L v ,d = Y[( vI + D L y ,d,i) 
i=l 

Tz,Gn,L y ,d — Vz : G D ;Ly,dHz,Gn,Ly,d , Tz,G D ,L y ,d = Vz : Go,L v ,dHz,Gn,L y ,d (7.4) 

where K is again the diagonal matrix with diagonal element w l where I is the number of vertices in the 
9 = 1 state for the corresponding basis, and K is the diagonal matrix with diagonal element w l where t is 
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the number of vertices, excluding the (L y + l)-th vertex, in the q = 1 state. Notice that Vz,G D ,L y ,d does not 
include the factor (vl + Al for the single external vertex. We have 

Ly + l 

Z {Gd, L y x L x , q, v,w)= c {d} Tr[Tz^G D ,L y ATz,G D x y ,dy n ^} (7-5) 

Here only one Tz,G D .L y ,d is needed as the single external vertex in the q = 1 state should only be considered 
once. Compare matrices Tz } G D ,L y .d and Tz,G D ,L y ,di certain columns differ by a factor of w, corresponding to 
the partitions with the (L y + l)-th vertex in the q = 1 state. It is clear that the number of these columns is 
nzh(L y , d), i.e. the cross case in category (b) in the proof of Theorem lll.H It follows that there are nzh(L y ,d) 
eigenvalues of Tz,G D ,L y ,d equal to the corresponding eigenvalues of Tz,G D L y ,d multiplied by w, and the rest 
of the eigenvalues are the same. Let us denote the eigenvalues that are common to Tz,G D ,L y .d and Tz,G D .L y .d 
as usual as Xz.G D .L y .d,j, and the eigenvalues of Tz.G D ,L y .d that lack a factor of to as Xz,G D ,L y .d,j- 

Furthermore, the common eigenvalues Xz,G D .L y ,d,j at level d also appear as \z,G D ,L y ,d+l,j at level d + 1 
with < d < L y . This again is due to the fact that Vz,g d ,l v a does not include the factor (vl + DL y ,dx y +i), 
so that it does not matter if the (L y + l)-th vertex is assigned a color (connected to a black circle) or not. 
Denote the number of these common eigenvalues as nzh(Go, L y , d) with 1 < d < L y + 1. They are given by 

n Z h(G D ,L y ,d) = n Z h(L y ,d- 1) + n Z h(L v , d) (7.6) 

That is, they can be classified as being in either category (a) plus category (c) in the proof of Theorem lII.ll or 
the circle case in category (b) plus category (d) with d replaced by d — 1 in that proof. In Table [V] we list the 
first few numbers nzhiGn, L y , d) and their total sums Nzh,G D ,L y - In particular, the numbers nzh{Go, L yi 1) 
is 

n Zh (G D ,L y ,l) = ]T (^V* (7-7) 



nzh{GD, L y , L y ) = 3L y , nzhiGo, L y , L y + 1) — 1 (7-8) 
If we denote the generating function of nzh(Gn, L y ,l) as 

oo 

Ai(x)= n zh (G D , Ly, l)x L » +1 = x + 3x 2 + 10x 3 + 36x 4 + ... (7.9) 

L y =0 

then the generating function of nzh{GD, L y ,2), denoted as ^(x), is given by the convolution [34| 

oo 

A 2 (x) = Y n zh {G D ,L y , 2)x L y +1 = [A^x)] 2 = x 2 + 6x 3 + 29x 4 + 132a; 5 + ... (7.10) 

L y =l 

In general, the generating function of nzh(Go, L y , d) is given as 

oo 

A d (x)= n Zh (G D ,L y ,d)x L * +1 = [A 1 (x)] d (7.11) 

L y —d— 1 

From Eq. (|7.6p . it is clear that 

N zh ,G D ,L y = 2N zh .L y - nzh(L y ,0) (7.12) 

It follows that 

Nzh t Go,L y +i = 5N Z h,Go,L y - n Z h{G D ,L y , 1) (7-13) 
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TABLE V: Table of numbers nzh{L y ,Gr>,d) and their sums, Nzh,G D ,L y ■ Blank entries are zero. 



L y \d 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Nzh,G D ,L v 


1 


3 


1 


















4 


2 


10 


6 


1 
















17 


3 


36 


29 


9 


1 














75 


4 


137 


132 


57 


12 


1 












339 


5 


543 


590 


315 


94 


15 


1 










1558 


6 


2219 


2628 


1629 


612 


140 


18 


1 








7247 


7 


9285 


11732 


8127 


3605 


1050 


195 


21 


1 






34016 


8 


39587 


52608 


39718 


19992 


6950 


1656 


259 


24 


1 




160795 


9 


171369 


237129 


191754 


106644 


42498 


12177 


2457 


332 


27 


1 


764388 



similar to Eq. (|2.23p . and Nzh,G D ,L y can be expressed as 



The multiplicity for the common eigenvalues with 1 < d is defined as 

d-1 



(7.14) 



k^(q) = K^(q-l) = c^+^ d - 1 ^Y,(- i y( 2d 1 3 )(l-^ (7-15) 

3=0 \ 3 / 

where k W (q) is given in Eq. (|7.3p . The first few of these coefficients are 
K™=q-1, = (q-l)(q-3) , 

= (q-l)(q-2)(q-4) , = (q - l)(q - 3)(q 2 - 6q + 7) (7.16) 

Collecting the common eigenvalues, we find that the partition function in Eq. (|7.5[) can be rewritten as 

Ly+l n Z h(L y ,GD,d) 

Z(G D ,L y xL x ,q,v,w) = £ {\z,G D ,L y ,dj) m 

d=i j=i 

L y n zh (L y ,d) 

+ w E £(d) E c^G D xy, dJ r (7.i7) 

d=0 3 = 1 

where nzh{L y ,d) is given in Theorem III. 11 



A. Determinants 



We find 



det(T z , GD r d ) = ^yyz^Ly+iA) 



n zh (L y A) 



io(l + J )(! + «) 



- L y n zh (L y ,d) 



(7.18) 
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for the full range < d < L y + 1. Comparing Vz, sq ,L y +i,d in Eq. (|2.27|) and Vz,G D ,L v ,d m Eq. (|7.4[) . one 
sees that Vz,G D ,i <j does not include the factor (vl + DL y .d,L y +i) for the self-dual square lattice, so that the 
power of v and (1 + q/v) has the factor L y rather than L y + 1. Comparing Hz,G D ,L y ,d and Hz,G D ,L y ,d in Eq. 
(|7.4p , one sees that the power of w for Hz,c D ,L y ,d has the factor rather than L y + 1 in Hz,G D -L y ,d- 

Taking into account that the multiplicity of each Xz,G D ,L y ,d,j is c^, it follows that the total determinant is 

det(T ZtGDtLy )= f[ [det{T z ,G D ,L y ,d)f d) = [v L »f v+1 U L " +1 (l + ^"(1 + v) L A " (7.19) 

d=0 L J 

B. Eigenvalue for d — L y + 1, A = Gd 

For d = L y + 1, each factor of Vz,G D ,£ s ,<i in Eq. (|7.4[) reduces to a scalar u, and Hz,G D ,L y ,d and Hz,G D .L y ,d 
become one. Both of the transfer matrices Tz,G D -L y ,d and Tz.G D .L y ,d for d = L y + 1 reduce to a scalar, namely 

^Z,G D ,L y .L y + l = V L « . 

C. Transfer Matrix for d = L y , A = G D 

The transfer matrix Tz,G D ,h y .L y has dimension 3-Lj, + 2. We obtain the following general formulas: 

{T Z ,Go,L y ,L y )lS = , {T Z ,G D ,L y ,L y )l,1 = V^W , {T Z ,G D ,L V ,L J2L B + 1,2L„ + 1 = V^ -1 ^ + <? - 1) (7.20) 



(T ZiGDiLy , Ly ) 2 j-i,2j-i =V L "- 1 (3v + q-l) for L y > 2 and 2 < j < L y 
{Tz,G D ,L v ,L y hj,2j = v^^wil + v) for 2 < j < L y + 1 

(7 1 Z,G D ,L H: L B )2L B +2 : 2L !/ + l = U^ 8 " 1 ^ + <? - 1) 

(7 1 z,g £ ,,l !; ,lJ2 J ,2 J -i =« i ^ 1 (2« + <z-l) for L y > 2 and 2 < j < L y 

(Tz,Go,L y ,L y hj-l,23 = for 2 < j < L y + 1 

(7V,G Dj L y ,i !/ )2j-l,2j+l = (7 1 Z,Gi,,L ! ,,L H -l)2j,2j + l = W H for 2 < j < L y 
(Tz,Gn,Ly,L y )2j+l,2j-l = (Tz,Gn,L y ,L y -l)2j+2,2j-l = V v for 1 < j < L y 

(Tz,Gn,L y ,L y ) j , j =v L y(l + v) for 2L y + 3 < j < ZL y + 2 

{Tz,G D ,L y ,L y )3,2L y +3 = (Tz,G D ,L y ,L y -l)4,,2L v +3 
{Tz,G D X y ,L y )2j-1.2L y +2+j = (Tz,G D ,L y ,L y -l)2j,2L y +2+j 

(Tz,G D ,L y ,L y )2j+l,2L v +2+j = {Tz,G D ,L y ,L y -l)2j+2,2L y +2+j = V Ly ~ l (l + v) for 2 < j < L y 



7.21) 
7.22) 
7.23) 
7.24) 
7.25) 
7.26) 
7.27) 
7.28) 



(7.29) 
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(Tz,G D ,L v ,L y )2L y +2+j,2j-l — {Tz,G D ,L y ,L y -l)2L y +2+j,2j+l 



for 1 < j < L y 



(7.30) 



with all other elements equal to zero. Tz,G D ,L y ,L y is the same as Tz,G D ,L y ,L y except for one element 

(T Z: Go,L y ,L y )2,2 = V L y. 

We illustrate these general formulas for the cases L y = 1 and L y = 2: 



T; 



Z,G D ,1,1 



/ v \ 

wv 

V X2 w V\ 

v X\ WVi Vi 

\ v 2 v 2 vv\ ) 



(7.31) 



G D ,2,2 



/ V 




















\ 





wv 

























Z3 


w 


V 





Vl 


Vl 







X2 


WV\ 


i) 





Vl 


Vl 








V 





X 2 


w 





Vl 








V 





Xi 


WVi 





Vl 


v 2 





V 2 











VVi 





V o 





V 2 





V 2 








VVl J 



(7.32) 



Thus, in general, the matrix Tz t G D -L y .L y is almost the same as Tz lS q,L y +i,L y except for the first two rows. It 
is obvious that Tz,g d ,l v .l v always has one eigenvalue equal to v Ly and one equal to v Ly w, and Tz.g d ,l v ,l v 
always has one eigenvalue equal to v Ly with multiplicity two. 

As corollaries of our general result for Tz,G D ,L y ,L y we calculate the determinant and trace: 



det{T z ,G D ,L y ,L y ) 
det(f z ,G D ,L y ,L y ) 



( v L y \3L y +2 
f y L y \3L y +2 



1 



tfl(l + J )(!+«) 



(7.33) 



which is the d — L y special case of (|7.18|) . and 



Tr(T z ,G D ,L y ,L y ) = v^- 1 L y {q - 1 + v 2 + Av) + w{v + L y + vL y ) 



Tr(T z ,G D ,L y ,L y ) 



,.Lv-l 



L y (q ~ 1 + v 2 + Av) + v + wL y (l + v) 



(7.34) 



VIII. SOME ILLUSTRATIVE CALCULATIONS 



A. Ly = 1 

We give some explicit results of the transfer matrices for cyclic strips in this section, beginning with the case 
L y — 1. We have exhibited the transfer matrix Tz, S q,isi m Eq. (|4.6p . and the quantity 1^,39,1,1 = v has been 
given above. The Potts model partition function Z(sq, L y x m, cyc.,q, v, w) was calculated for the the circuit 
graph in (43|. The results for eigenvalues agree. (The actual transfer matrices themselves are basis-dependent, 
and the basis used in (43j was different from ours, so the matrices are different, but the only part of the transfer 
matrices that enters into the partition function is the (powers of the) eigenvalues.) 
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B. Square-Lattice Strip, L y — 2 



The partition function for this case is given by the L y 
for d = is 



2 special case of Eq. (|2.3p . The transfer matrix 



T 



Z,sq,2,Q 



/(9- 


-1) 2 + 


3(9- 


l)v - 


f Zv 2 


X\W 


X\W 


W 2 V\ 


X 2 Vi 


\ 




-1) 2 4 


-2(?- 


l)t> 


+ v 2 


X\WV\ 


(q - l)w 


9 9 

w vf 


XiVi 






-1) 2 4 


-2(g- 


1> 


+ v 2 


(q - l)w 


XiVlW 


9 9 

w vf 


XlVl 






(9-1 


2 + (9 


-i; 


V 


(q - l)wvi 


(q - l)wvt 


w 2 vf 


(9 - l)«i 




V 




V 3 















V 2 Vl 


) 



(8.1) 



and the matrices Tz, S g,2,i & n d Tz.sq.2.1 have been given above. The result for the cyclic case Z(sq,2 x 
m, eye, q,v,w) agrees with Ref. 4J]. The partition function for the Mobius case follows from our general 
formulas also given above. 



C. Triangular-Lattice Strip, L y — 2 

We illustrate our results for the L y = 2 cyclic strip of the triangular lattice. We obtain 



f (q— l) 2 + 4(<7 - l)v + hv 2 + v 3 xiw 
(q - l) 2 + 3(q - l)v + 3v 2 + v 3 



Tz,, 



X1WV1 



X2W 

X\W 



W 2 Vi 



1 X\V\ 

(q — l) 2 + (q - l)v (q— l)wvf (q—l)wvi w 2 vf (q — l)vi 



(q — l) 2 + 2(q — l)v + v 2 (q — l)wvi xiV\w w 2 v 3 



(8.2) 



V 



V 2 V3 







wv 







v 2 v 1 v 2 ) 



The matrices Tz t tH,2,\ an d Tz t tn.2,2 were given above. 



D. Honeycomb-lattice Strip, L y = 2 



For the L y = 2 cyclic strip of the honeycomb lattice we calculate 



T, 



Z,hc,2,Q — 



I 




x\ 






w 


X2 






(9- 


- 1)31 




(9 - l)w 


W 2 V\ 


Xl 






(9- 


- 1)^1 


(9 - l)w 




W 2 Vi 


Xi 






(9 


-I) 2 


(g - l)w«i 


(q - l)wvi 


9 9 

w vf 


9-1 




\ 
















v 2 


) 



T; 



Z,sq,2,0 



(8.3) 



where Tz. sq ,2.o is given in Eq. 



The other matrices relevant for this strip were given above. 
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